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ABSTRACT 

The unusual pre-main-sequence binary star named KH 15D undergoes remarkably deep and long- 
lasting periodic eclipses. Some clues about the reason for these eclipses have come from the observed 
evolution of the system's light curve over the last century. Here we present UBVRI photometry 
of KH 15D based on photographic plates from various observatories, ranging in time from 1954 to 
1997. The system has been variable at the ^1 mag level since at least 1965. There is no evidence for 
color variations, with a typical limit of A(i3 — V)< 0.2 mag. We confirm some previously published 
results that were based on a smaller sample of plates: from approximately 1965 to 1990, the total flux 
was modulated with the 48-day orbital period of the binary, but the maximum flux was larger, the 
fractional variations were smaller, and the phase of minimum flux was shifted by almost a half-cycle 
relative to the modern light curve. All these results are consistent with the recently proposed theory 
that KH 15D is being occulted by an inclined, precessing, circumbinary ring. 

Subject headings: stars: pre-main sequence — stars: individual (KH 15D) — circumstellar matter — 
techniques: photometric 



1. INTRODUCTION 

The T Tauri star known as KH 15D (jK earns fc HerbstI 
Il998j) called attention to itself through a remarkable pat- 
tern of photometric variations. E very 48 days, the sys- 
tem dims by over 3 magnitudes ({Hamilton et alJl2001t 
IHerbst et al]l200l . It remains in this faint state for a 
duration that has grown from about 16 days in 1997, 
when the periodicity of the variations was first noticed, 
to over 24 days in 2004. Before 1999, the faint state was 
characterized by a short interruption of about 5 days dur- 
ing which the system returned to its bright state. Since 
then, these "rebrightening" events have subsided to near- 
u ndetectability. 

iHamilton et'all l|200lD and lHerbst et alJ l|2002D argued 
persuasively that these dramatic and unusual variations 
are due to eclipses by circumstellar dust that surrounds 
the visible K7 star, although the arrangement of the dust 
and the cause of the periodicity were unknown. At the 
time, it was not even clear whether the system was a 
single star or a multiple star system, and hence whether 
the relative velocity of the dusty structure and the visible 
star was due to the orbit al motion of the dust, or of the 
star ijHerbst et al.ll200^ . The mystery of the eclipses 
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of KH 15D motivated an intensiv e monitor ing cam- 
paign by observers aro und the world ijHerbst e t al. 200^ 
iBarsunova et al.*2004') , attracted other observers seeking 
clues in deep imaging ( Toku naga ct al. 2004) and spec- 
troscopy ifAgol et al . 2004; D eming et al.lll2 004'l. and also 
attracted t heorists hoping to model the circumstellar en- 
vironment lIBarge & Vitonll200l lG rinin fc Tambovtseval 
l2tffll . 

An important step in the clarification of the basic 
system properties was the analysis of archival photo- 
graphic plates of NGC 2264, the 3 Myr old star clus- 
ter in which KH 15D resides. Photographs from years 
prior to 1951 showed no evidence of eclipses deeper than 
1 mag (Winn et al. 20 0^. More interes tingly, photom- 
etry from 1967-1982 Ij Johnson fc WinnI 12004. hereafter 
Paper I) showed that eclipses were occurring with the 
same 48-day period as today, but with three striking 
differences: the system was brighter at all phases, the 
fractional variations in the total light were smaller, and 
the phase of minimum light was shifted relative to the 
modern light curve. 

Motivated by these findings, iWinn "eTaLl ll200l pro- 
posed that KH 15D is a binary star that is gradually 
being occulted by an opaque foreground screen, which is 
probably the edge of an inclined, precessing, circumbi- 
nary disk. C hiang fc Murrav-Clavl l|20041 had the same 
idea independently, and explored the dynamical issues 
further, arguing that the "disk" was likel y to be a fairly 
narrow ring. A radial velocity study by iJohnson et alJ 
|200?) confirmed that KH 15D is a binary star, and to 
the extent that the orbital parameters of the binary could 
be constrained by the radial velocities, the orbital pa- 
rameters agreed well with the predictions of I Winn et alJ 
( 2004). These developments have made KH 15D less 
mysterious than it once was, and they bring closer to re- 
ality the prospect of using the system to learn about disk 
evolution, the circumstellar environment of young stars. 
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and possibly even planet formation. 

A prerequisite for a complete model of the binary or- 
bit and the surrounding disk or ring is a light curve with 
reasonably complete time and phase coverage over the 
entire era of the periodic variations. Towards that end, 
we present additional photometry of KH 15D based on 
an analysis of photographic plates from several archives 
around the world, ranging in time from 1954 to 1997, 
with most of the data coming from years prior to 1985. 
Our plate selection and digitization process are described 
in § 121 and the photometric measurements and uncer- 
tainties are described in § 01 In § 01 we present the 
light curves and describe the evolution that have been 
observed over the past 50 years. In §|Slwe discuss these 
findings in the context of our model of the system. 

2. THE PLATES 

In Paper I, we presented results from an analysis of 52 
photographic plates in the archive of the Astrophysical 
Observatory of Asiago, in northern Italy. All of those 
plates were derived from observations with a single tele- 
scope, the 67/92 cm Schmidt telescope, and almost all of 
them were exposed with a red filter (RG5, or equivalently 
RG665) and a red-sensitive emulsion (I-N). We deliber- 
ately selected this subset of plates for the first phase of 
our archival study of KH 15D, because they formed a 
fairly dense time series with uniform characteristics, and 
because the effective band pass was approximately the 
same as the Cousins / band, in which most of the modern 
CCD-based photometry has been obtained. In addition, 
the photometry was simplified for the I-N/RG5 plates 
because the glare from a nearby B star is much smaller 
on red-sensitive plates than it is on the more commonly 
used blue-sensitive plates. 

Now that our photometric procedure has been tested 
thoroughly, we have expanded our analysis to include a 
heterogeneous collection of plates, obtained with differ- 
ent telescopes and effective band passes. The sample of 
87 plates that is considered in this paper is drawn from 
the following sources^: 

1. There are 26 plates from Asiago Observatory, the 
same archive from which the sample in Paper I 
was drawn (Barbicri ct al. 200^. But whereas 
the previously analyzed plates were taken with 
the 67/92 cm Schmidt telescope, the plates in 
the present sample were taken with the 40/50 cm 
Schmidt telescope and the 1.22 m Galilei telescope. 
As in Paper I, these plates were digitized at Asi- 
ago Observatory with a 1600 dpi, 14-bit commer- 
cial flatbed scanner. 

2. There are 35 plates that were taken with the 
105/150 cm Schmidt telescope at Kiso Observatory, 
in Japan, a facility of the University of Tokyo. The 
Kiso plates were digitized at the observatory with 
a PDS microdensitometer. 

3. There are 7 plates that were taken with the 
134/200 cm Schmidt telescope of the Thiiringer 
Landessternwarte Tautenburg, which is located 

® Most of these plates were identified by querying the Wide Field 
Plate Database, http://www.skyarchive.org 



near Jena, in Germany. They were digi- 
tized using the local high-precision plate scan- 
ner TPS at 2540 dpi and 12 bits pixer^. 
l|Brunzendorf fc MeusingeillT999ll . 

4. There are 7 plates that were taken with the 4 m 
Mayall telescope at Kitt Peak National Obser- 
vatory, which is located near Tucson, Arizona. 
The plates were shipped to the Harvard College 
Observatory for digitization with a 1600 dpi, 14 
bits pixel^^ commercial flatbed scanner. 

5. There are 5 plates taken with the 60 inch telescope 
of Mount Wilson Observatory, in the San Gabriel 
mountains of California, which were also shipped 
to the Harvard College Observatory for scanning. 

6. There are 4 plates that were taken with the 
2 m Ritchey-Chretien-Coude (RCC) telescope at 
the National Astronomical Observatory in Rozhen, 
Bulgaria. These plates were scanned with a 
1600 dpi, 14 bits pixel"^ scanner that is similar 
to the one used at Asiago Observatory. 

7. Finally, there are 3 plates from the Palomar Obser- 
vatory Sky Survey obtained with the 122/183 cm 
Oschin Schmidt telescope. These were scanned us- 
ing the USNO PMM plate scanner. The digitized 
images were downloaded from the Space Telescope 
Science Institute web site^". 

3. PHOTOMETRY 

The scanning process produces a digital record of the 
density of opaque grains in the photographic emulsion as 
a function of the position on the plate. As in Paper I, 
we used a form of profile photometry (also known as 
point-spread-function, or PSF, photometry) in order to 
estimate the apparent magnitude of KH 15D on each 
plate, relative to a grid of local reference stars. This 
procedure is described fully in Paper I and is summarized 
below. 

The relationship between the density enhancement, 
d{x, y), and the flux received from the star, S, is nonlin- 
ear, and for bright stars the density saturates at a maxi- 
mum value. Thus, an accurate PSF model must include a 
nonlinear resp onse. We employed the model described by 
iStetsonI l|1979() . in which the PSF model is a magnitude- 
dependent transformation of a two-dimensional Gaussian 
distribution, 

ds{x,y)^[dG{x,y)-^ + d;^,y'^\ (1) 

where dg^t and q describe the saturation properties of 
the plate, and dcix^y) is a Gaussian function of the co- 
ordinates X and y. The Stetson function has the desir- 
able limits ds — > dc for do ^ rfsat, and ds — > dsat for 
da > dsat- 

The model is fitted to the stellar image using a non- 
linear least-squares algorithm, and S is computed as the 
two-dimensional integral of dG{x,y). This flux scale is 
related to a standard magnitude via the expression 

m = Too - 2.5(1 + e)logio 5, (2) 

http://archive.stsci.edu/dss/ 
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Fig. 1. — BVRI photometry of KH 15D over the last 5 decades derived from archival photograhpic plates. 



where m is the apparent magnitude in the effective band 
pass produced by the emulsion and fiher combination, 
and e accounts for nonhnearity. 

For each plate, we selected the band pass in the 
Johnson-Morgan-Cousins UBVRI system that is the 
closest match to the emulsion/filter combination that was 
employed (see Table QJ. Then we estimated mo and e for 
that plate using the set of 53 reference stars that we iden- 
tified in Paper I. Those stars were selected by virtue of 
their proximity (they are all within a 20' box centered on 
KH 15D), their wide range of magnitudes and colors, and 
their relatively low level of variability. After measuring 
S for each of the reference stars, we found the values of 
mo and e that minimize the sum of the squared residuals 
between Eqn.[21and the magnitudes (in the app ropriate 
band pass) measured bv lFlaccomio et al.l l)1999|) . In Pa- 
per I, we noted that higher order nonlinear terms and 



color corrections were unhelpful in reducing the scatter 
in the magnitude relation. Wc found the same to be true 
for the present sample. 

For almost all the plates, these steps were carried out 
by an automatic procedure, but there were a few plates 
that required individual attention. Two of the Tauten- 
burg plates had an effective band pass approximating 
U band, but iFlaccomio et al.l 1)19991) did not measure U 
band magnitudes of our reference stars; for th ose plates, 
we used the U — B measurements of iPark et a l. (2000). 
On the Mt. Wilson plates, there was a strong gradient 
in the background near KH 15D, due to the glare of the 
nearby B star; for those plates, we found it necessary 
to account for this gradient by fitting a first-order two- 
dimensional polynomial surface to the background. Fi- 
nally, visual inspection of the POSSII IR plate revealed 
that KH 15D was much fainter than than any of our 
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Fig. 2. — Effective / band photometric time series of KH 15D, based on the assumption that the system has always had the same color 
indices as the modern uneclipsed state. 



reference stars. Rather than extrapolating to find the 
magnitude of KH 15D on this plate, we instead assign a 
firm lower limit of / > 15.5 based on the magnitude of 
the faintest of our reference stars. 

3.1. Estimation of Errors 

To estimate the uncertainty of each photometric mea- 
surement, we followed the same procedure as in Paper I. 
For each subset of plates obtained at the same telescope 
and with similar emulsion/filter combinations, we plot- 
ted the standard deviation in the resulting magnitudes 
of each star versus the mean magnitude of that star. The 
lower envelope of that relationship was taken to be the 
limiting uncertainty of our fitting procedure as a func- 
tion of magnitude, oi^rnin 

(/). We also took into account 
the relative quality of each plate by computing Cm, the 
scatter in the fit to Eqn. El and then we estimated the 
total uncertainty as 



UnC. = cr/,min(/) X 



(3) 



where crm,min 
that subset. 



is the minimum cr„j among all plates in 

4. RESULTS 

4.1. BVRI light curves 

The magnitudes of KH 15D are given in Table HJ The 
first four columns give the Julian Date, the observatory 
(using the designations defined in Table ^1, the clos- 
est corresponding UBVRI band pass, and the apparent 
magnitude with the the corresponding estimated uncer- 
tainty. The light curves for each band pass arc shown 
in Figure 1. For completeness, we have also plotted the 
measurements from Paper I. 



The BVRI light curves confirm two of the findings of 
Paper I: in the past, the system was formerly brighter 
and the fractional variations in total flux were smaller. 
The horizontal dashed lines in Figure 1 show the modern 
uneclipsed magnitudes of KH 15D. The light curves show 
that before about 1990, the system was often as much as 

1 mag brighter than it appears today. 

4.2. Test for chromatic variations 

iHamilton et"an l)200H) and lHerbst et all (|2002f) found 
that the flux variations that are observed today are 
nearly achromatic. Outside of eclipses, the V — I and 
V — R color indices do not vary by more than about 
0.05 mag. Even during the 3 mag fading events, the 
system hardly changes color, although at minimum light 
it does appear to be systematically bluer by wO.l mag 
in V — R. It would be interesting to know whether 
the flux variations in previous decades were also nearly 
achromatic, or whether they were accompanied by color 
changes. 

The measurements reported in Paper I were almost 
all in / band. We considered only a few plates with 
different band passes, and only one pair of plates with 
different band passes {B and R) that were taken nearly 
simultaneously (on 1974 December 15). The result was 
A(S - i?) = -0.14 ± 0.30, where A(B - R) is the dif- 
ference between the B — R index measured on that date, 
and the B — R index of the modern uneclipsed state. 

With the multi-band sample presented in this paper, 
we were able to test for color changes more extensively. 
There are 18 pairs of observations that took place within 

2 days of each other and that have different effective band 
passes. For each of those pairs. Table |21 lists the mean 
Julian date and the measured color index, C . The third 
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Fig. 3. — Phased light curves for different time periods. The filled circles are data based on the photographic plates described in this 
paper. The open circles (in the bottom panels) are modern CCD data, for which the error bars are smaller than the plot symbols. 



column lists AC, which the difference between C and 
the corresponding color index of the modern uneclipsed 
state^^. The final column gives the photometric phase, 
which is defined in § 14.41 

The result of this analysis is that there were no de- 
tectable color changes: all of the values of AC are con- 
sistent with zero, regardless of the photometric phase. 
Thus it appears as though the optical color indices of 
KH 15D have remained constant to within wO.2 mag over 
the past few decades. 

Here, for consistency, we have used the IFlaccomio et all 
11999) values for the une clipsed color indices. However, foUow- 
ine- IHamilton et alJ EHmt) . we have assumed that the B — V entry 
for KH 15D (star no. 391) in the catalog of Flaccomio et all II1999I) 
contains a typo, and should be 1.32 rather than 0.32. This makes 
the entry consistent with other published measurements (e.g. Park 
et al. 2000, who report B — V = 1.308) and with the color expected 
of a star of the observed K7 spectral class. 



4.3. An "effective" I band time series 

Under the assumption that the color indices have al- 
ways been exactly constant, we can put all the measure- 
ments on the same magnitude scale. The last column of 
Table |2] gives the "effective" /-band magnitude, which 
was computed under the assumptions U — I — 4.02, 
B - I = 2.94, V - I ^ 1.62, and i? - / = 0.80. This 
effective / band time series is shown in Figure [5] 



4.4. Evolution of the light curve 

We have divided the effective / band time series into 
different time intervals, bounded by the dotted lines in 
Figure |21 in order to examine how the time series has 
evolved over the past 5 decades. Since both the modern 
variations and the variations observed from 1967-1982 
(Paper I) are nearly periodic, it is useful to plot the mag- 
nitude as a function of the photometric phase 0, which 
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is defined as 

, (J.D. - 2,434,800) mod F 

^-p , (4) 

where P = 48.38 days is the orbital period that was de- 
termined bv l.Tohnson et a L ( 2004). We note that this def- 
inition of (j) differs from the photometric phase that was 
used in Paper I, because of the updated period. In Fig- 
ure 3 we have plotted phased light curves for each of the 
6 time intervals covered by our data, along with phased 
light curves from CCD-based data that were kindly pro- 
vided by W. Herbst and C. Hamilton. 

The phased light curves clarify the nature of the photo- 
metric variations. From about 1965 to 1991, the system 
alternated between a bright state and a faint state with 
a 48-day period, just as the system does today. How- 
ever, the modern bright state is approximately 1 mag 
fainter than the pre-1991 bright state. In fact, the mod- 
ern bright state has approximately the same brightness 
as the pre-1991 faint state. The other striking feature 
of the phased light curves is that (j)niim the phase corre- 
sponding to the midpoint of the faint state, has changed 
by nearly 0.5. Formerly it occurred at 4'min ~ 0.32. By 
1995-96, it had shifted to (j>min — 0.8, although it is worth 
remembering that in that year and the few years to fol- 
low, the "faint state" included a short-lived rebrightening 
event. 

These findings confirm and extend the results of Pa- 
per I, which were based almost entirely on the 1965-73 
data. The enlarged sample indicates that the shallow, 
phase-shifted eclipses were occurring until at least 1991. 
It also shows that the variations were smaller and less 
regular in the years before 1965. Between 1958 and 
1965, there is some evidence for a faint state centered 
at 0min ~ 0.7, but the fractional variation appears to be 
0.5 mag, and the transition between the bright and 
faint states is more gradual. Before 1958, no particular 
pattern of variations is evident. 

5. DISCUSSION 

IWinn et a,].l ll20fil and IChiang Mnrrav-Cl^ ll200l 

hypothesized that KH 15D is a pair of stars in a mutual 
orbit with high eccentricity, and the "eclipses" are actu- 
ally occultations by a circumbinary disk. A full anal- 
ysis of this model should include all the available in- 
formation, including the modern photometry, the pho- 
tographic pliotom£tr^jand__the radial velocity measure- 
ments bv l.Tohnson et al.l l|2004f) . This will be the subject 
of a future paper. Here we simply describe how the basic 
properties of the light curves presented in this paper can 
be explained by the model, and point out one aspect of 
the new data that has not yet been explained. 

In the model, we are viewing an eccentric binary sys- 
tem nearly edge-on. Today we can see only one of the 
two stars; we refer to the visible star as A, and the in- 
visible one as B. The reason B is invisible is that its en- 
tire orbit is in the shadow of a circumbinary disk. Also 
shadowed is a portion of A's orbit, near periapse. The 
3 mag eclipses occur when A passes into the shadow, 
once per period. The reason for the decadal evolution 
of the light curve is that the line of nodes of the disk is 
regressing, in response to the gravitational torque pro- 
duced by the stars. This causes the shadow of the disk 



to sweep across the system. As a progressively larger 
portion of A's orbit is shadowed, the eclipses increase in 
duration. Extrapolating backwards in time, star A used 
to be completely exposed, and only a portion of star B's 
orbit (near apoapse) was shadowed. In that era, we ob- 
served periodic occultations of star B. 

Thus, we interpret the historical light curve as follows. 
The periodic alternation between the bright and faint 
state is due to the periodic occultation of star B. The 
total flux was formerly larger, and the fractional vari- 
ations were formerly smaller, because the steady light 
of star A was visible at all orbital phases. Finally, the 
reason for the phase shift is that the previous eclipses 
occurred around apoapse, whereas the modern eclipses 
occur around periapse. 

Extrapolating even further backwards in time, one 
would expect to reach an era in which neither star was 
shadowed at any phase, and there were no photometric 
variations due to occultations. The new data presented 
in this paper show a somewhat more complicated picture. 
From 1951 to 1965, the system varied at the «0.5 mag 
level. This is certainly smaller than the 1 mag periodic 
variations observed between 1965 and 1991, but it is not 
zero. It is not clear whether these early variations were 
periodic, or erratic. The 1958-1965 data do seem to obey 
the 48-day periodicity, but the data are too sparse and 
the variations are too small to be certain. 

In closing, we wish to note that although we have at- 
tempted to find as many useful photographic plates of 
NGC 2264 as possible, we do not claim that our collec- 
tion is complete. Undoubtedly there are additional plates 
in other observatory archives. Additional photometry 
would be useful for a quantitative determination of the 
orbital parameters of the binary star and its putative 
circumbinary disk, especially if those data were taken 
during the years 1980-1995 where the time coverage of 
our light curves is not very complete. 
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TABLE 1 

The plate sample 



Observatory 


Name of 


Range of 


Number of 






Approximate 


Scanning 


(Designation) 


Telescope 


Years 


Plates 


Filter 


Emulsion 


Band 


Method 




ZLn/'^n rm "^rhmirlt 




20 


Rr;5 ROfifi^ 


I-N 


J 


I. ItxL UCU 


Aqiao-n rAmi9n^ 


±.Zi^ 111 Vjtdlllt,! 


1955 


g 




103a-O 




FlathnH 

J. itxL LJl^(_l 




J-UiJ/ ±0\J L.1I1 xVloU OL.llllilU.lj 








Ila-O 




PDS 










Y 










4 


r;r:4Q5 


lla-D 


Y 










1 




103a-F 


R 










2 




Ila-F 


R 










4 






R 










g 




I-N 


J 




Kitt Peak (KPNO) 


4.0 m Mayall 


1978-1981 


2 


RG610 


IIIa-F 


R 


Flatbed 






2 


RnfiQ5 


IV-N 


J 












00385 


Illa-J 


Y 










1 


GG495 


Ila-O 


B 










1 


GG385 


Illa-J 


B 




Mt. Wilson (MTW) 


60 in 


1954 


3 


GGll 


103a-D 


V 


Flatbed 






2 


GG13 


103a-O 


B 




Palomar (POSS) 


122/183 cm Oschin Schmidt 


1951-1996 


1 


RP2444 


103a-E 


R 


PMM 




1 


RG610 


IIIa-F 


R 










1 


RG9 


IV-N 


I 




Rozhen NAO (ROZ) 


2.0 m RCC 


1990 


4 


GG385 


ZU21 


B 


Flatbed 


Tautenburg (TAU) 


134/200 cm Schmidt 


1972-1997 


4 


GG13 


ZU 


V 


TPS 








2 


UG2 


ZU 


U 










1 


RGl 


103a-E 


R 





TABLE 2 

Photometric Measurements of KH 15D 



Julian Date 


Observatory 


Band Pass 


m [mag] 


/eff [mag]'' 


2433658.9 


POSS 


R 


14.2 ±0.4 


13.4 


2434801.0 


MTW 


V 


15.4 ±0.2 


13.7 


2434802.0 


MTW 


V 


15.2 ± 0.2 


13.5 


2434805.0 


MTW 


V 


15.5 ± 0.1 


13.9 


2434806.0 


MTW 


B 


16.8 ± 0.4 


13.9 


2434833.0 


MTW 


B 


16.7 ±0.4 


13.7 


2435452.5 


ASI120 


B 


16.35 ± 0.08 


13.41 


2435452.6 


ASI120 


B 


16.32 ± 0.07 


13.38 


2435458.5 


ASI120 


B 


16.43 ± 0.08 


13.49 


2435458.6 


ASI120 


B 


16.29 ± 0.08 


13.35 


2435459.5 


ASI120 


B 


16.32 ± 0.06 


13.38 


2435460.4 


ASI120 


B 


16.31 ±0.07 


13.37 


2436846.6 


ASI040 


I 


13.4 ±0.2 


13.4 


2436849.7 


ASI040 


I 


13.5 ±0.10 


13.5 


2436876.5 


ASI040 


I 


13.4 ±0.1 


13.4 


2436881.6 


ASI040 


I 


13.5 ±0.1 


13.5 


2436935.4 


ASI040 


I 


13.49 ± 0.09 


13.49 


2436952.4 


ASI040 


I 


13.6 ±0.1 


13.6 


2436954.4 


ASI040 


I 


13.7 ±0.1 


13.7 


2437259.5 


ASI040 


I 


13.6 ± 0.2 


13.6 


2437348.3 


ASI040 


I 


13.75 ± 0.08 


13.75 


2437349.3 


ASI040 


I 


14.0 ±0.1 


14.0 


2437350.3 


ASI040 


I 


13.8 ± 0.1 


13.8 


2437369.3 


ASI040 


I 


13.8 ± 0.1 


13.8 


2437373.3 


ASI040 


I 


13.26 ± 0.07 


13.26 


2437379.3 


ASI040 


I 


13.4 ±0.1 


13.4 


2437562.6 


ASI040 


I 


13.8 ±0.10 


13.8 


2437584.6 


ASI040 


I 


13.44 ±0.09 


13.44 


2437619.6 


ASI040 


I 


13.5 ± 0.1 


13.5 


2437649.5 


ASI040 


I 


13.9 ±0.1 


13.9 


2437659.6 


ASI040 


I 


14.1 ±0.2 


14.1 


2440578.4 


ASI040 


I 


13.17 ±0.07 


13.17 


2441335.0 


TAU 


R 


15.2 ±0.1 


14.4 


2442726.3 


KISO 


B 


16.1 ±0.2 


13.2 


2442747.0 


KISO 


I 


14.6 ±0.10 


14.6 


2442747.1 


KISO 


R 


15.2 ±0.1 


14.4 


2442747.2 


KISO 


R 


15.5 ±0.1 


14.7 


2442811.1 


KISO 


V 


16.0 ±0.2 


14.3 


2442812.1 


KISO 


B 


17.4 ±0.3 


14.5 


2442859.9 


KISO 


V 


15.7 ±0.2 


14.1 


2442868.9 


KISO 


V 


15.0 ±0.1 


13.4 


2442870.9 


KISO 


B 


16.5 ± 0.2 


13.6 


2442871.0 


KISO 


V 


15.1 ± 0.1 


13.5 


2443073.3 


KISO 


B 


16.7 ±0.3 


13.7 


2443085.2 


KISO 


V 


16.1 ±0.2 


14.5 


2443085.3 


KISO 


V 


15.9 ± 0.3 


14.3 


2443085.3 


KISO 


R 


15.2 ± 0.1 


14.4 


2443113.2 


KISO 


R 


14.29 ± 0.08 


13.49 


2443113.3 


KISO 


R 


14.25 ± 0.08 


13.45 


2443192.9 


KISO 


I 


14.48 ± 0.09 


14.48 


2443210.0 


KISO 


I 


13.62 ± 0.05 


13.62 


2443211.0 


KISO 


I 


13.63 ± 0.06 


13.63 


2443213.0 


KISO 


B 


16.5 ±0.2 


13.5 


2443215.0 


KISO 


B 


16.6 ±0.3 


13.7 


2443217.0 


KISO 


B 


16.7 ±0.3 


13.7 


2443217.0 


KISO 


R 


14.3 ± 0.1 


13.5 


2443241.0 


KISO 


B 


17.5 ± 0.4 


14.5 


2443244.0 


KISO 


B 


17.2 ± 0.4 


14.2 


2443244.0 


KISO 


V 


15.8 ± 0.2 


14.1 


2443769.3 


KISO 


I 


14.42 ± 0.08 


14.42 


2443821.0 


KPNO 


R 


15.25 ± 0.09 


14.45 


2443821.1 


KISO 


V 


16.0 ± 0.1 


14.4 


2443833.2 


KISO 


I 


13.66 ± 0.06 


13.66 


2443848.0 


KPNO 


I 


13.87 ±0.09 


13.87 


2443949.0 


KPNO 


B 


17.4 ±0.1 


14.4 


2444644.0 


KPNO 


V 


16.0 ±0.2 


14.4 


2444644.0 


KPNO 


R 


15.19 ±0.08 


14.39 


2444644.0 


KPNO 


I 


14.43 ± 0.07 


14.43 


2444644.0 


KPNO 


B 


17.4 ±0.1 


14.4 


2444942.6 


KISO 


V 


15.9 ±0.2 


14.2 


2445996.7 


KISO 


V 


16.2 ±0.2 


14.6 


2446003.6 


KISO 


B 


17.4 ±0.6 


14.5 


2446494.4 


KISO 


B 


16.7 ±0.2 


13.7 


2446496.4 


KISO 


V 


15.1 ±0.1 


13.5 


2446776.2 


KISO 


R 


15.6 ± 0.2 


14.8 


2447592.0 


TAU 


B 


17.4 ±0.2 


14.5 


2447592.0 


TAU 


B 


17.3 ± 0.2 


14.4 


2447837.9 


POSS 


R 


15.4 ±0.2 


14.6 


2447947.0 


ROZ 


B 


17.01 ± 0.08 


14.07 


0/(/17Q/17 n 




R 


1« «Q 4- n (17 


1 Q Qt; 
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TABLE 3 
Color Measurements of KH15D 



J.D. 




Index 




C 




AC'' 


0b 


2434801, 


.5 


R 

rj — 


V 


1. 


,3 ± 0. 


,2 


0.0 


0.03 


2442747, 


.0 


-ft. — 


T 
1 


0. 


.6 ± 0. 


1 


-0.1 


0.26 


2442747, 


.0 


D 

JX — 


T 
1 


0. 


,9 ± 0. 


,2 


0.1 


0.26 


2442811, 


.0 


B — 


V 


1. 


.5 ± 0. 


,2 


0.2 


0.58 


2442870, 


.0 


D 

LJ — 


V 


1. 


.5 ± 0. 


1 


0.2 


0.80 


2442870, 


.0 


B - 


V 


1. 


.4±0. 


1 


0.1 


0.80 


2442871, 


.0 


B - 


V 


1. 


,6 ± 0. 


1 


0.2 


0.83 


2443085, 


.2 


V - 


R 


0. 


.9 ± 0. 


,2 


0.1 


0.25 


2443085, 


.2 


V - 


R 


0. 


,7±0. 


,2 


-0.1 


0.25 


2443212, 


.0 


B - 


I 


2. 


,8 ± 0. 


1 


-0.1 


0.87 


2443217, 


.0 


B - 


R 


2. 


.4±0. 


,2 


0.2 


0.98 


2443244, 


.0 


B - 


V 


1. 


.4±0. 


,2 


0.1 


0.53 


2443821, 


.0 


V - 


R 


0. 


,9 ± 0. 


1 


0.1 


0.46 


2444644, 


.0 


B - 


V 


1. 


.4±0. 


,2 


0.1 


0.47 


2444644, 


.0 


V - 


R 


0. 


.8 ± 0. 


,2 


0.0 


0.47 


2444644, 


.0 


R- 


I 


0. 


,8 ± 0. 


1 


0.0 


0.47 


2446495, 


.5 


B - 


V 


1. 


,5 ± 0. 


1 


0.2 


0.74 


2450481, 


.5 




B 


1. 


.3 ± 0. 


,4 


0.2 


0.13 



^AC = Cmcasur cd ~ C'o wherc Co is the color 
index reported bv lFlaccomio eT al. ( 1999 ) (or, in 
the ca se of the single U — B entry, nPark et alj 
|2003)). 

^Binary orbital phase, as defined in Eqn.|l] 



